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PREFACE. 



The purpose of this book is to present those portions 
of Differential Equations which are most used by engineers. 
It was written with reference to the needs of two classes of 
students: those who take the subject as a preparation for 
the study of Mechanics, and students in Electrical Engineer- 
ing who are already well grounded in the physical ideas in- 
volved in the electrical problems. The former take the first 
three chapters, omitting the electrical problems; the latter 
complete the book. 

I am indebted to Edwards' Integral Calculus, and still 
more to the text-books of W. W. Johnson and D. A. Mur- 
ray. The symbols and equations of the electrical problems 
are practically the same as in Bedell and Crehore's Alternat- 
ing Currents. In my own classes the books will be used to- 
gether. 

The examples were prepared with the double purpose 
of illustrating the text, and of affording a review of some of 
the more important opertions of the Calculus. 

Answers are given and stress is laid on the importance 
of verifying results. 

In the problems, the physical meaning and use of the 
arbitrary constants are kept prominent. 

I am under obligations to Dr. H. W. Kuhn for sugges- 
tions and corrections. 
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DIFFERENTIAL EQUATIONS. 



CHAPTER I. 

DIFFERENTIAL EQUATIONS DEFINED. 

1. A differential equation between two or more variables i& 
any equation containing differentials of the variables. 

dy = fix J dx. dy = (x^—a') dx, 

U{y)dy = f{,\ , ax . fdy = x'^dx, 

fi(^,y)dy = fix)dx, ix' + y'')dy = logxdx, 

are some examples of differential equations. 

2. The order of a differential equation is determined by that 

of the highest derivative in it, and its degree is the highest degree 

of its derivatives. 

d^y ^ dy _^ 



dx'^ dx 
is an equation of the second order and the first degree, while 

\dx ) ^ dx 
is one of the first order and second degree. 

3. The ordinary integral calculus deals with the differential 
equation, in which the derivative is expressed explicitly as a func- 
tion of one of the variables, which may be regarded as the inde- 
pendent variable. In the equations (/y=(jr* — a^)dx, dy=f(x)dx, 
X may be taken as the independent variable, and we have in the first 



dy _ 



dx 
and, in the second 



Kx'^-a\ 



^^ --f^Y 



dx 
In each case the derivative is a function of x only. In like manner 

the equation 

dy = f(y)dx 
may be written 

dx _ 1 

dy f{y)' 



2 INTEGRATION CONSTANTS. 

in which case the derivative with respect to y as the independent 
variable is a function of y only. 
In the same way^ the equation 

^-■^^ 

is a problem of ordinary calculus, and may be solved by integrating 
twice with x as the independent variable. 
On the other hand, the equation 

can not be integrated directly, since y can not be taken as the in- 
dependent variable. 

INTEGRATION CONSTANTS. 

4. Since the derivative of a constant is zero, a constant may 
vanish each time an expression is differentiated. Hence, when per- 
forming the reverse operation, a constant must be introduced on 
integrating an expression of the first order; two constants, on 
integrating an expression of the second order, and so forth. 

We really introduce a constant in the ordinary integral calculus, 
when we integrate between definite limits, but it vanishes in the 
process of subtracting the value of the expression at the lower limit 
from that at the upper. 

The value of a constant of integration depends upon the condi- 
tions of the problem in which it occurs. In a general solution, it 
should be of such form that it mav have any value from H-ooto — oo. 
Generally, it is expressed as an algebraic quantity as C, or, if the 
result is more symmetrical as O or O, etc. Sometimes it is more 
convenient to give it some other form, as log C, tan C, tan— *C, 
etc., which fulfill the condition that they may have any value from 
-footo — 00, Such an expression as sin C can not be used unless 
it is known that the conditions of the problem are such that the 
constant must lie between + 1 and — 1, or unless C is imaginary. 

In the differential equation 

xdx'-\-ydy^^^, 
the solution is evidently 



X 



2 1/2 



f-2:L=c' 



2 ' 2 

which may be written more symmetrically 

Considered geometrically, this is the equation of a circle of radius C 
Consider the differential equation of the motion of a body hav- 
ing a uniform acceleration g. This may be written 

dv 
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from which 

The constant C, in this case, is the value of the velocity when 
* = 0, in other words, it is the initial velocity. 

A constant of integration, or, as it is generally called, an ar- 
bitrary constant, has a definite value in a particular problem, but 
may have an entirely different value in anpther similar one. 

The integral of a differential equation with the number of arbi- 
trary constants corresponding to the order of the equation is called 
the general solution, complete integral, or complete primitive. 

EQUATIONS OF THE FIRST ORDER AND FIRST DEGREE. 

VARIABLES SEPARABLE. 

5. The simplest form of a differential equation is one in which 

the variabfes may be separated, so as to bring it into the form 

f (x) dx = f (y) dy, which may then be integrated directly. For 

example : 

ix + a)dy=i^+b)dx. 

This may be written 

y-\-b x+a 
The integral of which is 

log(y + &)=log(^ + a)+C'; 

or log(y + b)=z \og{x + a) + log C ; 

whence 

y^=CiX'{'a) — b. 

Examples. 
1. 4-=^' 3y* = 4x' + C. 

dx 

dv — 

^ 3. .-r2L=ydx, y=Ce '^ 

X -^ 

V 4. xydy={f + \){x-^l) dx, ^xHy»+l) = C^^. 

5. J^=sec ydXy sin y=^+C 

X 2 

6. -^=4L, t^ny=^+C. 

sec*jK dx 2 

v' 7. (1+/) dx=il + x') dy, tan-l>^=tan-l;ir+t*n-lC; 

otherwise, y=^Ll—^ 
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8. »fl—3^dx=^l^x^dy. 






>=;rVl— C«-|-CV1— ^*, 



otherwise, C^ =x^—y^—y^\ — xK 



'=^'<^y- 



FORMATION OF A DIFFERENTIAL EQUATION. 

6. A differential equation may be formed from any equation 
containing a single arbitrary constant by differentiating once and 
eliminating the constant between the result and the original equa- 
tion. If the equation contains two arbitrary constants, differentiate 
twice, and so forth. Take the answer of problem 2 of the pre- 
ceding list of examples 

x'y + Cy + B = 0, (1) 

Differentiating once we get 

Sx^y dx+(x^ + C) dy = 0. (2) 

Eliminating C between equations (1) and (2) we get 



dy ^^,y_ 



dx 



(3) 



which is a differential equation free from the arbitrary constant 
C. If the constant stands alone, it vanishes on differentiation, so 
that it is not necessary to eliminate. For instance, divide equation 
(1) above by y, and we have 

y 

which differentiated gives 

Zx^dx—^dy^O, 

a differential equation free from the arbitrary constant C. 

Examples. 

Check the results of the examples of Art. 5 by the method 
given above. 

HOMOGENEOUS EQUATIONS. 

7. Any equation of the first order and first degree, which is 
homogenous in x and y may be put into the separable form by sub- 
stituting vx=^y, or vy=^x, where z/ is a variable. Take the equa- 
tion 

(:i^-f-/) dx-k-xydy — ^ (1) 

Let y = Vjr, hence dy=^v dx-^ x dv. 
Substituting in (1) 

.r*(l-|-2^) dx -\- ^viv dx '\' x dv) (2) 

Dividing by ^r'and collecting 

(l+2z/')(/^-f t/jr(fv = 0. (3) 
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This may be separated into 

_^+^=0. (4) 

Log AT-h-i-logd -h2z/8) =log C. ( 5 ) 

4 

jr*(l + 2v*)=C* (6) 

Substituting z;=-:L.we get 

jr* + 2;ry = C*. (7) 

The general homogenous equation is of the form 
{Ax^ -^Bx^-^y-^Cx'^-^y^ .. . .-\-/^y'*)dx 

+ (^'Ar« -f-^'^«-'j/ -h Cat" -*y +^>« )rfy=0. 

Substituting y = z/^, rfy = vdx -\-xdv, we get 

+xn {A'-\-B'v+ C'lP' -\-R'v^ ) {vdX'\-xdv) =0. 

Dividing by x^ and collecting 

[A-\-{A'+B)v+{B'-\-C)7^ -^R'vn+i]dx 

^{A'-\-B'v+C'v^ +R'vn)xdv=0. 

d,x . ■ A'-{-B'vJ^C'vt -\-R'vn ^^^^ 

X '^A-\-(A'-\-B^v-{-{B'+C)v^ -\-R'v»+i 

From this it is evident that any homogenous equation is rendered 
separable by substituting y = vx, and, from the symmetry of the 
terms, the same is true for x=:vy. 

Examples. 
1. (y'^2xy+2x')dy=y'dx, ' y—x=C(y^—2xy). 

V 2. {x+iijy+xy)dy=ydx, Cy=J[^ ^-^. 



3. {y-{-^x^-\-y^)dx=xdy, Cx^-y-\-^x^-\-y\ 

Otherwise, Cjt' — 2Cy = l. 

4. Check the answers of the above problems by differentiating 
and eliminating the constant of integration. 

NON- HOMOGENEOUS EQUATIONS OF THE FIRST DEGREE IN X AND Y. 

8. Equations of the first degree in x and y which are not 
homogenous are of the form 

{ax^hy + c)dx^- (a'^r+fry + cOrfy^O. (1) 

Such equations may generally be transformed to the form 
(ax' + hy)dx' + (aV+ h'y')dy '= C, (2) 

which is homogenous and may be solved by the methods of the 
preceding article. In equation O) let 

x=^ x' -\-h, y = y' -f k, hence djt — dx\ and dy =. dy' , (3) 
where h and k are undetermined constants. 



} 
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Equation (1) then becomes 
iax'+ by'+ ah + bk + c) dx'+ 

{a'x'+ hy+ a'h + h'k + c')dy' = 0. (4) 
If now we assign to h and k such values that 

ah-\-hk + c = 

: (5) 

a'h + b'k + c'= 
Equation (4) becomes (2). 

It is necessary, then, to solve the homogenous equation (2), 
which, neglecting the prime, is the same as (1) with the con- 
stants omitted, and substitute from (3) x' =^x — h, y' =y — k. 
Take the equation 

(x + 2y + l)dx-(x'^l)dy=:0, (1) 

Substituting and dropping the primes 

ix+2y)dx — xdy = 0, (2) 

which solved as a homogenous equation gives 

Cx' = y + x. (3) 

For the constants h and k we have 

h + 2k + l=0, 

(4) 

A— 1 = 0, 

Whence h = l, and k = — 1. 

Substituting x — 1 for jr and y + 1 for y in equation (3) we get 

C(x-iy = y + x. 

Examples. 

1. (2x + y + l)dx+(x + 2y^4)dy=0, 

x' + xy-\-y^ + x — 4y = C, 

2. ix + 2y + S)dx=i2x + y + S)dy, 

C(x--yy = x + y + 2, 
^ 3. (4X'-y + 6)dx+(x — 4y + 5)dy = 0, 

ix^y + 2y(x + yr=C, 

4. (x--y)dx+(x + l)dy = 0, ^^^^^^..^^ 

9. If in a differential equation of this kind, — r-=— r;- the 

method given above fails, for the equations for finding h and k 

are of the form 

ah + bk + c = 0, 

rah + rbk + c' = 0, 

the loci of which represent parallel straight lines when c' is not 

equal to re, and the same straight line when c' = re. In the first 

case, the values of h and k are infinite, and in the second case, they 
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are indeterminate. These are easily solved, however, by a dif- 
ferent method. For the indeterminate case divide by (ajr+6y+^). 
Take the equation 

(2:r — 4y + 6)fl?:r— (3Ar — 63; + 9)dy = 0. 
Dividing and transposing, we have 

2dx = My, 
Hence 2;r = 3y + C 

For the first case, let ax + hy^^z^ a new variable. The method 
is shown by the problem below: 

{x + y + 2)dx+{2x-\-2y + \)dy = (^, (1) 

Let ^ + y = ^, hence dx = ds — dy. 

i2 + 2)(d2 — dy) + (2s + l)dy = 0, (2) 

(s + 2)d2+(s — l)dy = 0, (3) 

£r + 3; + 3 1og(^-l)=C, (4) 

x + 2y + S\og(x + y-l)=C. (5) 

Examples. 

1. (x + y)dx—{x + y + l) dy = 0, 

2x — 2y + C = log(2x + 2y + l), 

2. (2x + y + S)dx+(6x + Sy + 9)dy = 0, x + Sy = C. 

3. ix + 2y + S)dx— i2x + 4y + l)dy = 0, 

^ — 2y + log(4^ + 8y + 7)* =C 

Miscellaneous Examples, 

1. (x+l)dy=(y + 2)dx, Cix+l)=y + 2. 

x — l 

2. {X'\-y)dy={y + \) dx, C+Cy=eyTi 

4. ix + y)dx+i2x + 2y + l)dy = 0, 

logCix + y + l) = 2y + x. 

EXACT DIFFERENTIAL EQUATIONS. 

10. An exact differential equation is one which is formed by dif- 

ierentiation alone without elimination of constants, or reduction. 

X dy-\- y dx =^0 

is an exact differential equation, being formed by differentiating 

xy = C. 
If we differentiate 

nx,y)=C, 
we get 
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dx ' dy 

as the standard form of an exact differential equation. This may 

be written 

Mdx + Ndy = 0, 

where M is the partial derivative of the function with respect to 

jr and N is the partial derivative with respect to y. 

An equation of the form 

Mdx + Ndy = ^, 
is exact if 

dM _ dN 

dy dx ' 
For 

M=-^/ix,yh N=-^/{x,.yi, 

that is: --^ — is obtained by differentiating f{x,y) with respect 
dy 

dN 

t© X and then with respect to y ; while -^ — is obtained from 

ax 

the same function by differentiating with respect to y and then 
with respect to x, and as it makes no difference in what order 
these operations are performed, the expressions are equal. To il- 
lustrate, take the equation 

x^y + Sxy = C, 
differentiating 

(2xy + 3y)dx + (x' + 3x)dy = 0, 

an exact differential equation. 

Differentiating (2xy-\-^y) with respect to y we get 

^{2xy-{-3y)=2x+3=^. 
dy dy 

Differentiating (x'"^ + '^x) with respect to x, 

-^{x^+S.r)=2x-tS=^. 
dx dx 

SOLUTION OF AN EXACT DIFFERENTIAL EQUATION. 

11. Integrate M dx, regarding j as a constant and :r as a var- 
iable, then integrate the terms of A^ dy which do not contain :r, 
regarding 3; as a variable, and add the results. 

For example we will solve the differential equation 
Ox'' + ^xy)dx + {2x' + 2y)dy = 0. 
This is an exact differential equation for 

_^(3.r24-4;f r)=4;r, and -^(2x^+2y)=^. 
dy ' dx 
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Integrating (3:r*+4jry)rfjr, regarding x as variable, we get 

The only term of {^Ix^ -\-^y)dy which does not contain x is 2ydy, 
the integral of which is y. 
Adding, we have for our solution 

jr» + 2^3' + y = G 
which may be verified by differentiation. 

The reason for this method of integration is evident from 
the following considerations: Any function of x and y is made 
up of three se;ts of terms: (a) those containing x only, (b) 
those containing y only, and (c) those containing both x and y. 
On differentiating to form an exact differential equation (a) pro- 
duces terms of Mdx which, of course, do not contain y, (b) pro- 
duces terms of Ndy which do not contain x, while (c) produces 
terms of Mdx which contain y and terms of Ndy which con- 
tain X.. If we integrate the terms of Mdx which do not contain 
y we get the terms which contain x only. If we integrate the. 
terms of Ndy which do not contain x we get the terms which con- 
tain y only, while if we integrate either the terms of Mdx which 
do contain y or the terms oiNdy which contain x we get the terms 
in X and y. 

Examples. 

Show that each of the following differential equations is exact 
and solve. 

1. (2xy' + 2y)dx+ (3^^ + 2x + 2y)dy = 0, 

x'y' + 2xy + y' = C. 

2. 2x dx + exy dy + Sy^dx = 0, - x^ + Sxy^ = C. 

y y^ 2 y 2 
xdx-\-ydy dx 

5. yxy-^dx-^xy\o%xdy=^, x^ = C. 

EQUATIONS WHICH MAY BE RENDERED EXACT. 

12. Sometimes an equation may be made exact by multiply- 
ing by a factor. The equation 

x^dx -i-xy dy = 

is not exact, but if it is multiplied by — it becomes 

X 

xdx-\-ydy = Oy 
which is exact. The factor which renders an equation exact is 
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called an integrating factor. Integrating factors may frequently 
be found by inspection. 

There are a number of methods for finding integrating factors. 
We will give one which is applicable to cases where the factor is 
a function of x only or of y only. [For other methods, see Mur- 
ray's Differential Equations.] Let the equation be 

and let R be the integrating factor. The equation 

RW dx-\-RN' dy^^, 
is exact, and 

^ {RM')=^-^{RN'), (1) 



dy dxr 

Suppose R is SL function of x only. Since 

^^0. and 4^=^^ 
dy ox^ dx 

^ bM'_ RbN , N' dR .g. 

by dx dx 

R n\ dy dx J ' ' 

If R is, as we have assumed, a function of x only, the coefficient 
of dx in (3) is a function of x and is directly integrable. Take 
the differential equation 

{x'f — y)dx+ (x'y + x)dy = 0, 
Assuming that R is 3. function of x we have from (3) 

dR __ (2x^y—l—Ax:^y—\)dx __ 2dx 
R x^y-\-x X 

Hence, R— — - is the integrating factor. 



x^ 



If the expression 



1 f dM' _ dN' ^ 



N'\ dy dx 

is not a function of x^ we know that the assumption that i? is a 
function of x only, is incorrect. We may then assume that R is 
a function of y, and we get from (1) 

As the expression on the parenthesis in equations (3) and (4) is 
the same, the work of finding the integrating factor, whe^ it is a 
function of one variable only, is not great. Take as another ex- 
ample. 

(Zx^y^ —— ") ^■^+ {x^y'\-l)dy=Oy 
dy dx X 
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If we attempt tp divide the term to the right of the equality sign 
by the coefficient of dy, we find that the result is not a function 
of X only, consequently the integrating factor is not a function of 
X, The term is divisible by the coefficient of dx, the quotient being 

1 

y 

Hence from (4) 

dR _ dy „„, 1 
—--=——^1-, and — 

R y y 

is the integrating factor. The exact differential equation is 

(zx^y—l-Sdx-^-lx^^ 1^ dy =0. 

Examples. 

Find integrating factors, solve the equations and verify the 
results. 
/ 1. y(l+y)rf^— (2;r + 2;ry-y)rfy = 0, 

>' 1 + Cy 
2. (3;»ry — 3y + 3:r+l)rf;ir+ U— l)Vy=:0, 

-^^ x—\ {x—i)^ 
^' 3. ^x'dx + xdy — ^ydx = 0. y = C^' — 3:r*. 

LINEAR EQUATIONS. 

13. A linear differential equation is one in which the dependent 
variable and its derivatives are of the first degree. The standard 
form of linear differential equation of the first order where y is the 
dependent variable and x the independent variable, is 

dy,^Py=Q^ 



dx 
P and Q being functions of x only. 

This equation may always be made exact. To prove this 
statement and find the integrating factor, arrange the terms and 
multiply by R, which is the assumed integrating factor. 

R{Py — Q)dx^Rdy = 0. (1) 

This equation being exact 

^ R{Py-Q)=4^ (2) 



dy dx 

If /? is a function of x only (2) becomes 

dx 
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SPtix 

Substituting this value of R in (1) and rearranging the terms in 
the order of the standard form, we have 

e V^'^dy + Pye% ^^""dx = Qe^ ^^""dx, ( 3 ) 

The derivative with respect to x of the coefficient of dy is 

Pe^ ^, and the derivative with respect to y of the coefficient of 

^^ .« »..«>, - - , consequently equation (3) is exact. 

To solve (3) integrate the coefficient of dy, regarding y as 
the variable, and those terms of the coefficient of dx, which do 
not contain y, regarding x as the variable. The final result is 

yeV^'^^QeVdx^C. , (4) 

The differential of the first term of (4) gives the first two terms 
of (3), while the differential of the two terms to the right of the 
equality sign of (4) is the last term of (3). Take as an example, 

4-^^^^\ (1) 

ax X 
jPdx =C ^^ =log X, 

Equation (3) becomes xdy •\-ydx^ix^dx, 
the solution of which is 

x^ ^ x^ C 

xy=^^C, or>^=^+ii-. 

4 4 X 

When y is an explicit function of x and of the first degree, 

it is frequently convenient to verify the result by substituting this 

value of y in the original equation. For instance, if we substitue 

x^ C 

in equation (1) above the value — + — =JV» we get 



dy 



dx X 4 x^ 4: x^ 



Examples. 



1. ll2L+2xy=4x, y=2—Ce-^. 

dx 

3. ^— jK=a, y=Ce^—a. 

dx 
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4. ^+j/=^, y^ x—\-\-Ce-x ^ 
ax 

5. ^^Zy^x, y=\^^\,J^Ce-^ 
ax 3 9 

6. ^-\-y=8in Xf y=JL~{sinx-—cosx)+Ce—^^ 
ax 2 

7. ^+ay=sm bx, a sin ^^-^ cos bx_^^^^ 
dx a*+^ 

To integrate e<^^ sin &;»: rf;»r, we make use of the method of 
integration by parts. Since d(uv) =udv-\-vdu 

j udv =uv — ( vdu. 
If we let ^^ dx = dv and sin hx = u, this becomes 

/e^^ sin bxdx = — e^^ sin dx — — ie^*^ cos bxdx, ( 1 ) 
a a J 

We may operate in the same way on j e^^ cos bxdx and get as the 

final term — —5- | ^«^ sin bxdx. As this is the same integral as the 

one with which we began, it may be transposed to the left of the 
equality sign and added. The solution will then be obtained by 

2 AS 

dividing by^i-i — . Another method is to integrate the original ex- 
a^ 

pression by parts, letting sin bxdx—dv and e^x-=u. This gives 

j e<^x sin bxdx^^J^ax cos \ix +-^ Je*^ cos bxdx, ( 2 > 

Multiplying (1) by a and (2) by h and adding 

( a2+ ^ J e<^^ sin bxdx=e<^^ ( a sin bx-^-b cos ^^ ). 

To put the answer into the second form requires a trigonometric 
transformation ; a sin bx — h cos bx may be put into the form of the 
sine of the difference of two angles by dividing by a factor r such 

tliat JL, equals the cosine, and — . equals the sine of the same angle. 
r r 

To find the value of r, 

a sin ^;r — bcosbx_^/ a . b , \ 1 
o , Lo 1 ^ sm bx — „ cos bx I „ , . 

from which the second form of the answer is obtained. 
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g, Siiij^-^!L-|-4'-lc<»jr=— ;r«, Cos ^=:£l-^rjr. 

ydx 






10, ^-1— iL__^^2 j'tj-l)=Z1^4l-fr. 

11. dy'\-—^:==fix^xdx, yxx^\\ —x*)= ^ 



EQUATIONS WHICH MAY BE TfcAXSFOHifED INTO LIXEA&S. 

14. An equation of the form 



dx 

i» a linear. If we divide by y^ it becomes 

which is erjuivalent to 

-^ ^(y^-n)-^Pyl-^=Q. 

1 — n dx 

Multiplying by d — n) this is 

dx 
a linear in which y^ — « replaces the y of the standard form. 
As an example 

y ^xy—\xy'' 



dx 

dx 

d{y~^) -1 

"11^ --J—xy =—4x, 

dx 

_ X2 X^ 



y e i =4e '^ +C. 
(4-f Ce- -^ )y—l=0. 



Examples. 

1. _^l'_+>'ta.i;.-l5Ili:, J'»=l+<rcos«;r. 

a.v y 

dv .»* -I ;i— 1 3 ^ ^+1 
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dx y ' 5 

Miscellaneous Examples. 
Solve and verify results : 

^ 2. xdy + 2ydx = x*dx, 6y = x* + Cx-'^. 

S. (x + y)dx-{- (2x + 2y — S)dy = 0, 

x + 2y + S\og(x + y'-3) r=C. 

4. 2y(/;r— (Ar=^ + 2^y — y)(fy=0, yC^ + y) = C(ar — y). 

5. fdx + xYdy + idx = 0, y = —+ 2 tan - ^-^-\-C. 

>y 6. B dy -\- y tan X dx = y^cot3in X dx, 

y^sec ;ir = log tan^ + C. 



7. xdx^ i^x^-\-y/^dX'\-ydy — 0, x+ i^x^+y^ =C. 

8. -^+e^yz=e^y y = l-\-Ce-^ 

dx 

9. (iy + x^)dx = xdy, y = Cx*^x^. 
10. ^Tan y + — ^—^dx -\- x s^d^ydy = 0, x taxi y+tan- 1 x—C. 

PROBLEMS. 

15. 1. Find the equation of the curve for which the tangent at 
any point is perpendicular to the radius vector. 

The slope of the tangent to any curve is -^ . The tangent of 

dx 

the angle which the radius vector makes with the X axis is-i-- 

X 

Hence the differential equation is 

dy X 

dx y * 

the solution of which is 

x' + f = C, 

a circle of radius C. The constant C may have any value from 

to 00, unless limited by conditions other than those stated above. 

The result means that in any circle, whatever the radius, the tangent 

is perpendicular to the radius. 

2. Find the curve for which the tangent of the angle which the 

radius vector makes with the x axis is twice the slope of the 

tangent to the curve. 

The parabolas y*=zCx, 
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3. Find the curve for which the tangent and the radius vector 
make supplementary angles with the x axis. 

Any rectangular hyperbola of the form xy = C. 

SOME DIRECT CURRENT CLBCTBICAL PROBLEM S. 

16. In these problems we will use the following symbols: 

C capacity, 

R resistance, 

L inductance, 

q quantity at any instant, 

i current at any instant, 

e potential difference at any instant, 

e base of natural logarithm, 

Q maximum quantity, 

/ maximum current, 

E maximum E.M.F. or potential difference. 

1. The terminals of a charged condenser are connected through 
a non-inductive resistance. Find the expression for the charge at 
any instant. 

The current strength is -^ aid e = -^. The current is also 

equal to the rate of discharge, that is, i —. — —i-, hence, we have the 

at 

differential equation 

9 — _ dQ i\\ 

RC di 

Equation (1) is the case of variables separable, and its solution is 



RC 

q= C^e '^'' . VI) 

The constant C\ is obtained from the condition that when ^ = 0, 

t_ 

RC 
q = (2, and consequently ^ = 1. Hence, Cj = Q, and (2) be- 

comes 



q = Qe ''''. (3> 

Equation (3) may be written 

log^ J- = — i_, or log^ _^=_J_. ( 4 > 

^ Q RC' ^ q RC ' 

A condenser of 2 micro-farads capacity is discharged through 

a resistance of 100,000 ohms. What portion of the charge remains 

in the condenser after .1 sec, .2 sec, .5 sec, 1 sec, and 2 seconds? 
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In this problem, ^C = 10« X 2 X 10-* = .2, and -JL = 5. The 
work may be conveniently arranged as shown below : 

Time, ^g^ -y '^Sio'y ^Ki(r^ -^ 

.1 .5 .21715 —1.78285 .6065 

.2 1.0 .43429 —1.56571 .3679' 

.5 2.5 1.08573 —2.91427 .0821 

1.0 5.0 2.17147 —3.82853 .0067 

2.0 10.0 4.34294 —5.65706 .0000 

Show that the charge in the condenser decreases in a geometricaP 
progression. 

If we can measure a quantity equal to one ten-thousandth of 
the original charge, how long will it take to discharge this con- 
denser completely as far as our measurements will enable us to 
determine? 1.842 seconds* 

2. A condenser is connected suddenly through a non-inductive 
resistance to a source of constant E.M.F.. Find the charge at any 
time t. 

If E is the E.M.F. of the source, the effective potential dif- 
ference at any instant is £ — e, hence 

E-e^Ri=R% (in 

at 

is the differential equation required. The solution of (2) is 

t^ 

q^CE=C^e ^^. (3> 

Substituting for CE its value Q, and inserting the condition that 
^ = when ^ = 0, to get the arbitrary constant Q, we find as the. 
final result 



= c(i-. ^'0. 



(4> 



3. In problems (1) and (2) find e and t. 

/_ 

From {\)e = Ee , 

Q L E i_ / 



Prom (2)e 
2 D B 



= ^(1-. '0. 
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What is the meaning of the negative sign before the first ex- 
pression for current? What assumpiion in the derivation of the 
formula determines this sign? 

4- The terminals of a constant p^xemial dynamo are suddenly 
connected through an inductive resistance. Find the expression for 
the current in terms of the time, the resistance of the entire circuit, 
and the inductance of the entire circuit. 

The E.M.F. must be equal to the drop of potential due to the 
resistance, together with the counter E.M.F. due to inductance. 
Hence we have the physical equation 

E^Ri + J^. (1) 

at 

Hence, dividing by L we get the standard linear form 

di Ri E .^^ 

~dr^~T = ~L' ^''' 

of which the solution is 

E Rt 

t = -^ -\- C^e z. , { 3 ) 

The constant C\ is determined from the fact that when / := 0, i ^), 

F 
hence, Cj = — — = — /. The final equation is then 

R 

Rt 



i = l(l-e 



5. The terminals of a coil of inductance L, in which a current 
/ is flowing, are suddenly joined together through a non-inductive 
resistance and, at the same instant, the battery supplj'ing the cur- 
rent is disconnected. Determine the current after an interval of t 

second. 

_^ 

t ^ le , 

R being the sum of the resistances of the circuit. 

If the inductance is 20 milH-henrys, the resistance of the coil 2 
ohms, the non-inductive resistance 98 ohms, and the initial potential 
difference 2 volts, find the current and potential difference after 
.0001 sec, .0005 sec, and .001 sec 

i .606 .0821 and .007 amp. 

e 59.44 8.04 and .66 volts. 

6. A cell having an E.M.F. of 2 volts and negligible internal 
resistance is connected to a cell of 1 ohm resistance and .01 henry 
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inductance. After the current has became steady 9 ohms non-induc- 
tive resistance is suddenly cut into the circuit. Derive an expres- 
sion for the current, and compute the current, the potential dif- 
ference of the terminals of the cell and of the terminals of th« 
non-inductive resistance after .0001 sec. 

^• = .2 + 1.8^r-l«»^ 
Potential difference of terminals of coil 16.29 volts. 

ALTERNATING CURRENT PROBLEMS. 

17. The typical form of an alternating E. M. F. is expressed 
« = £ sin w^, when w is equal to 2^ times the number of complete 
periods per second, and t is measured from a position in which the 
e is passing through zero in the positive direction. 

1. Circuits having resistance and inductance. 

Our differential equation after dividing by L becomes 

di , R ' E . . 

—+-^t=—s.n «/. (1, 

The solution of this linear equation is 

E 



ie ^ =-j-j e ^ sin oJt dt + C. ( 2 ) 

This equation is the same as Ex. 7 of Art. 13. Performing the 
integration in the same way we get 

Ri 
i — £(ig sin <at — L<a cos <at) . r-T'lL r^ , 

/e2+.L2a,2 + ^^ ' f3) 

which reduces to 



E 



inf w/ — fan -1-^)+^^ ^, (4) 



t = SI 

V i?' + L V 

Tan — 1-^ is the angle of lag, and "^ R^ -{- L V is the impedance. 
R 

C is the arbitrary constant. 

We will find the constant C for a numerical problem. Given 
L = .045 henry, i^ = 24 ohms, £^= 150 volts and w 400. From this 
we get the impedance to be 30. We can not make the assumption 
that t = when / = as in the preceding direct current problems, 
for we have already made the assumption that / — when e=0. 
The time of one complete period is .0157 sec. Curve I of Fig. 1 
represents E. M. F. for a little more than two and one-half peri- 
ods, the abscissas being time and the ordinates electromotive force. 
One division, representing .002 second, corresponds to 45°, 50' of arc. 
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Curve II of the same Fig. represents 



(sin«^ — 0). 



The angle of lag <t> is 36°-52', corresponding to .00161 second of 
time. This curve first crosses the zero line .00161 to the right of 
the corresponding point of the E. M. F. curve. 

The constant C depends upon the point where the circuit is 
closed. In the figure we have plotted the case when it is closed 

.003 second after initial time. We have (wf — 0) =3r-53', ^= 

1.6, and find 6V-1 • = —2.641, and C = — 2.641<?l-6 = —12.93. 

Rt 
Having C, we may now plot the curve Ce ^ . This is curve III 
of Fig. 1. 
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Figure 1. 



In reality we need only that part of the curve beyond t = .003, 

80 we might as well have nsed Ce ^ and plotted from A to the 

right. The current curve IV is the sum of curves II and III. At 

the instant of closing the circuit, as we have assumed in finding C, 

the current is zero and curve III is just as far below the zero line 

as curve II is above. The current curve IV then starts at zero at 

.003 sec. and lies below curve II until it reaches about .01 sec, where 

curve III becomes negligible. From that point it coincides with 

the sine curve II. The effect of the exponential term in this case is 

appreciable for about the time of one alternation. 

Rt 

For given values of R and L the curve Ce ^ has the same form 
whatever the value of C. For different times of closing the circuit 
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it is only necessary to displace curve III to the right or left so as 
to make it pass through the point whose abscissa is the time of 
closing and whose ordinate is that of curve II with the opposite sign. 
For instance, suppose that the circuit is closed when t is .006 sec. 
The ordinate of II is 4.9, and curve III must pass through the point 
whose co-ordinates are .006 and — 4.9. 

EXAMPLES. 

1. Find the value of C for the following cases : When the cir- 
cuit is closed at .006 second, .009 second, .012 second, and .015 second. 

First case, C = — 120 
Last case, C = 11920. 

24/ 

2. Construct the curve be 045 on a sheet of thin cross-sec- 
tion paper and construct curve II on a second sheet. Combine the 
two to find the current curves for the cases when the circuit is 

closed at .002 sec, .004 sec, .0096 sec, and .012 second. 

Rt^ 

The curve Ce ^ represents the fall of a direct current of ini- 
tial value C in a circuit of resistance R and inductance L. We may 
imagine that, at the time of closing the circuit, the alternating cur- 

E 
rent ./-— - . ^ sin (w^ — 0) begins at once at its full value, but that 
^IC -\- L(>r 

there exists with it a direct current in the opposite direction, which 

at that instant, is exactly equal to the alternating current, and which 

dies away as if the alternating current did not exist. 

Circuits having resistance and capacity. 

An alternating E. M. F. is suddenly applied to a circuit consist- 
ing of a condenser in series with a non-inductive resistance. Find 
the current at any time. 

The E. M. F. equation is 

RiJ^A^ = E sin «f , ( 1 ) 

Substituting for i its value-^ we get the linear 

The student will show that the solution of this linear reduces to 
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From (3) we get by differentiation 

E 



/^^3i:sinl(u,. + tan-i^)+C,. CR. (4) 



where C, = ^^ 



^C 



Tan — ^ _„ is the angle of advance, and may be represented by 
CR*a 

0. If the time of closing the circuit be represented by /, we have, 
since ^ = 0, 

If « == 400, R = 500 ohms, C = 10 microfarads, and e = 150 sin 
«^ construct q and i curves the circuit being closed when t .003 sec. 






CHAPTER II 

EQUATIONS OF THE FIRST ORDER AND 
HIGHER DEGREES. 

It is customary, when dealing with equations of higher degree, 

to express the derivative, -j^, by the single letter p. An equation 

ax 

of higher degree then, is a function of x, y, and p. There are three 

principal classes of these equations: 

1. Equations solvable for p. 

2. Equations solvable for y. 

3. Equations solvable for x. 

EQUATIONS SOLVABLE FOR P. 

19. When an equation can be solved for p without giving 
terms in x and y too compllicated for integration, this method should, 
in general, be used, as it gives us one or more equations of the first 
degree to which the methods of the preceding chapter apply, for 
example : 

(gy-3;rg.+2**=0. (1) 



Solving for p : 



ax ax 



The solutions of these equations are: 



x^ 



y = x^ + C, andy = ^ + C. (2) 

Each of these is a solution of the differential equation. Both may 
be combined to form a single solution : 

iy-C-x^) {y-C-:^)=0. (3) 

As this is the solution of a differential equation of the first order 
we need only one arbitrary constant in the result, and accordingly 
we drop the prime from the C in the second term. 

We will show now that (3) is the solution of (1) by perform- 
ing the operation of differentiating and eliminating the constant. 
Differentiating (3) and dividing by dx: 

(y—C—^){p—x)\{y—C—^) {p — 'lx) = 0. (4) 

^-C = ^iMzii£), (5) 

■^ 2(2/ — 3^) ' 

23 
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Substitmmg iht yalae of (y—C) in equatiDn i3i and 
wc get 

jr'igj- — /' ^/ — jri _^ 

whence we get ^^ — 2jr) (p — x) = 
which is equatron (Ij. 

EXAMPLES. 

Solve the following difiterential equations and Tcriiy the an- 
swers; 

EQUATIONS SOLVABLE FOR Y. 

20. An equation which can not be conveniently solved for p, 
tn2Ly often be solved for y, the factors being in the form of 

y^fU.P). (1) 

To integrate such an expression first differentiate it, and divide by dr. 
Since -i = /, we get 

a differential equation which does not contain y. Regarded as an 
equation in x and pA^) is of the first degree. Accordingly it may 
be solved by the methods of the preceding chapter, giving a result 
as a function of x, p, and C. The final result is obtained by elim- 
inating p between this result and the original equation (1). Take 
the equation 

2y = px + p\ * (1) 

Differentiating and dividing by dx 

2/-/+(;r + 3/2)^. (2) 

f-f = ^^' <^> 

a linear equation, the solution of which is 

x = ^p'+Cp. (4) 

we must now eliminate p between (4 J and (1). If we attempt to 
do this by substitution we shall get results involving surds. It is 
better therefore to use another method. If we multiply (1) by 3 
and (4) by /> and substract, we get 

C/>' — 4M+6y = 0. (5) 

We have eliminated />''and have now two equations (4) and (5) of 
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the second degree. We will eliminate />* between these last two 

equations and get 

p{C^Vlx)=Cx-\-\%y, (6) 

Substituting the value of p from (6) in equation (4) or (5) we get 

2Cf — CV + SeCxy — 16x^ + 108:y = 0. 

EXAMPLES. 

Solve and verify results. 



' 1. >' + -(fy=o. { 



/ 2. y^-\-/^ = 0, J,, (or 4- (7)2 + 4 = 0. 

(log^)2 + 21og/ = 2;ir + C 






(— 1 -f- V2^+ C-l-1) . 



EQUATIONS SOI<VABI*K FOR X. 

21. When an equation can not conveniently be solved for/ or^ 

it may sometimes be solved for Jt:. In that case after solving for x^ 

dx 1 
we diflPerentiate and divide by dy. Calling ^^-X"* ^^ ^^^® * ^^^ 

equation in y^ p^ and -^ . This is integrated and p is eliminated 

dy 

between the result and the original equation. 

Take the equation 

This may be written 

x-{-^ = 0. 

Differentiating, dividing by dy and solving the linear thus formed 

Substituting the value of /— ^ in (4) in equation (1) we get 

xy + 20^*+ r«j^ + 1 = 0. <^) 

That (5) is the integral of (1) may be shown in the usual way by 
differentiating once and eliminating C. 

1 dx 
In problems of this kind, instead of using — for ---to introduce 

t dv 
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another new variable, the letter q is generally nsed in works on dif- 

ferential equations for>-^ The problem is then exactly the same as 

dy 

that of Art 20, except that x and > change places. Using this nota- 

tion, equation ( 1 ) of the preceding paragraph becomes 

x+g»y = 0. (1> 

Differentiating and dividing by dy 

y + 3^y-r2«y»^ = 0. (2) 

ay 



From (3) we get the linear 



dy^7,y 2>«' 

3 



(3) 



? ~^ X ^ (4> 

qy =y -\-C. \^' 

Substituting (4) in (1) we get the same result as in the preceding 
case. 

EXAMPLES. 

1. ;r('^y+y=0. 4r4-(logC>)« = 0. 



2. 



claihaut's equation. 
22. An equation of the form 

y = px + fip). (1) 

is called Clairaut's form. Treating it as an equation solvable for 
y we get 

P=P + (x + f (P))-^, (2) 

(:r + /'(p))-g- = 0. (3) 

Equation (3) is satisfied by jr + /'(/>) = 0,and-^ = 0. 

From the latter we get p ^ C, which substituted in (1) gives as 

the general solution 

y = CX + f(C). (4) 

All that is necessary to be done in solving an equation of Clairaut's 
form is to substitute C for p in the original differential equation. 
There still remains the other factor of (3), 

^ + np) = 0, (5) 

which, regarded as a function of p, and x does not contain any dif- 
ferentials. A solution may be obtained therefore without any inte- 
gration by merely eliminating p between (5) and (1). The result 
thus obtained, since there is no integration, does not contain any 
arbitrary constants, and is called a singular solution. We have 
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then two solutions of Clairaut's equation; a singular solution with- 
out arbitrary constants and a general solution with an arbitrary 
constant. 
Take as an example 

y = xp+^. (1> 

Differentiating and dividing by dx 
hence 



C'+i) 



P\dp _, 



dx 



= 0. (3> 



From / = 0, we obtain the general solution, 
dx 



y=Cx + ^, (4> 

4 



If we eliminate p between x + -E^ = and equation (1), we get the 



^^^^p -. 

singular solution 

^ + 31 = 0. (5) 

If instead .of eliminating p we integrate the equation 

^ 2dx 

we get 

^ + y = C^, (6) 

which is the same as the singular solution with the addition of 
the constant Ci. 

The meaning of these results in a geometrical problem is seen 
from Fig. 2. The straight lines represent the equation 

y = Cx+—, 

for values of C differing by .2. For instance, the line 4 — 4 corre- 
sponds to the case where C is .4, the line 6 — 6, to the case where C 
is .6, etc. 

The singular solution, jr^ + y ^ 0, is the equation of the para- 

bola (not drawn in the figure) to which all lines y =^Cx -{- ——are 

4 

tangent; that is, the singular solution is the envelope of the lines 

representing the general solution. 

The equation x^ -{- y =^ C\, represent parabolas with vertices 

displaced upwards a distance C\. We will show that any one of 

these parabolas is the locus of the intersection of lines, 

y=Cx Jr ^. 
4 
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for which C differs by a definite amount. In the figure we have 
drawn the parabola for the intersection of the lines for which the 
values of C differ by unity. One point on this curve is the 
intersection of the line — for which C is with the line 10 — 10 
for which C is 1. Another point is the intersection of 2 — 2 with 
12 — 12. To find the equation of this curve let y = Cx "h_^and 

4 

y =■ iC -\- \)x -^ -i — i — L be the equations of the intersecting 
lines. Eliminating C we get x^ -\- y — 



16 



Show that the locus of the intersection of the straight lines 

y = Cx + — r- ^or which the values of C differ by a constant 
4 

a 

amount a is the parabola x^ -\- y =^ Ci, where Ci = -^. 

16 




Uiz 



Figure 2. 



MISCELLANEOUS EXAMPLES. 



1. y = px+p\ 

dy 



2. y=x^ 
ax 



S. x=y -^ + 



4. 



dy 
dy \2 



~d7' 



{-%-y^^-^'^^='' 



6. Tan-l-^ = ;tr-l, 
dx 



General solution, y = Cx -\- C; 
Singular solution, 4jr* + 27/ = 0. 

y = Cx — log C; 



and y •= \ -\- log x, 

or x^ — ^y = Cj. 

x^ 

y -\-x = Cdt— . 

y + log cos {x — 1 ) = C, 



CHAPTER III. 

DIFFERENTIAL EQUATIONS OF ORDER HIGHER 

THAN THE FIRST. 

23. We will take up the following forms of differential equa- 
tions of higher orders: 

1. Equations with y wanting, 

2. Equations with x wanting, 

3. Linear equations with constant coefficients. 

EQUATIONS WITH Y WANTING. 

24. These equations are of the form 

and are solved by simple integration, a constant being inserted 
each time an integration is performed. Take the equation 

d^ ~ 
Multiplying by dx this becomes 

Hence 



dx 4 

which integrated again gives 

25. This simple form of differential equation has many appli- 
cations. One of these is the case of a moving body. If .y is the 

ds 
distance of a body from some point taken as the origin, —-- is the 

dt 

velocity and , . is the acceleration. 

Problem. A body moves with a uniform acceleration a, find the 
velocity and position after / seconds. 

The constant Ci is the initial velocity, and d represents the initial 
displacement, that is, the position of the body with reference to the 
origin at zero time. 

29 
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Another important application is that of a beam with a trans- 
verse load. The differential equation of an elastic beam is 

EIJ^ = M. 

E being the modulus of elasticity of the material, /, the moment 
of inertia of the cross-section with respect to the neutral axis, and 
M the moment of the applied forces with respect to the neutral axis. 

A beam is fixed horizontally in a wall (Fig. 3), and projects 
a distance, /, to the right- A load, P, rests on the free end. Find 
the equation of the elastic curve and the deflection of the free end. 

A section distant x from the wall is / — x from the free end, 
and the moment at the section due to the load P is — P (/ — x). 

The negative sign is used because the slope of the beam— ^ becomes 
less, algebraically, as x increases. Hence the differential equation is 




Figure 3. 

To determine the arbitrary constant, Ci, we have the condition 

that the beam is horizontal at the wall. Hence — ^ = when jt = 0, 

dx 

which substituted in equation (2) gives Cx = 0. 

Integrating again, we get 

£/y = - P(J^ - ^^ + C, (3) 

and from the condition that at the wall where ;r ^ 0, ;y = 0, we get 
C2 = 0. The equation of the elastic curve is then 

At the end where .r = / the deflection is 

4 Pn 

-^ 3£"/ 

Solve this problem integrating (/ — x)dx as a single quantity. 
That is, write equation (2) 
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31 



EIJ^ = J^(/ — ^)2 + C. 
Show that the resuhs are the same, but the constants different. 




K\\\N 



't0fi 



Figure 4. 

Solve the same problem with the beam projecting to the left of 
the wall (Fig. 4), taking the origin at O, the position of the free 
end of the beam before it was deflected. 



EQUATIONS WITH X WANTING. 

26. Take first the equation of the second order 



d-y 



dx' 

calling — iL = /) as in Chapter II, this is 
ax 



(1) 



dp 



dx 



f(y), 



(2) 



To solve multiply dy dy, whence 

Pdp = f(y)dy. (3) 

This integrates to 

-^ = ff(y)dy-]^C. (4) 

Equation (4) is of the first order and second degree, and is integrated 
the second time by the methods of the preceding chapter. 



Examples. 



1. 



''^ = - 5y, 



dx^ 



2. -^>- = ^. 



X =--=sin-lI_--^-hC8. 
X = log {y + y/j,2 _|_ ^2 ) + (Tg. 



or^r = sinA — 1^-- + ^a- 

27. Example (1) has some interesting applications. One of these 
is the case of mass vibrating on a spring. When the distortion of 
a spring is inside the elastic limit, the force required to produce a 
given displacement is proportional to the displacement. If it 
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requires a fore, P, to produce a displacement of one unit, the force 
producing a displacement y is Py. If a body attached to a spring 
is displaced a distance y, the spring exerts on this body a force 
— Py^ the negative sign being used because the force is in a direc- 
tion opposite to the displacement. Since the product of the mass 
of a body, multiplied by its acceleration in any direction, measures 
the resultant force acting on the body in that direction we have 

In order that equation (1) may hold it is necessary of course that 
the units of mass and force be such that unit force produces unit 
acceleration in unit mass. If the mass is in grams, the force must 
he in dynes. If the mass is in pounds, the force must he in pound- 
als. If it is desired to use the pound force as the unit, the corre- 
sponding mass is the Engineer's Unit of 32.2 pounds, for which the 
name of gee-pound has heen suggested hy Mauer. 

ds 

Calling -— . V, and integrating as shown in Art. 26, we get 

_^ii = -^+Q. (2) 

To obtain the constant Ci we will call maximum displacement a. 

Pa^ 
Consequently when ;y = a, z^ is zero and Cx = — -— . Hence we have 

^ = ^(a'-y) (3) 

It may be easily shown that (3) is an energy equation. — s~^* 

the kinetic energy of the moving mass. — f-- is the work done in 

Pa^ 
displacing it a distance y, and -— - is the work in displacing to a 

distance, a, from the origin. Transposing the second term, we get 

Mi^ , Py^ __ Pa* 



2 2 2 ' 

which states that the sum of the kinetic and potential energies is a 
constant, and is equal to the work done in producing the maximum 
deflection. In fact, we might have begun by writing this equation 
from the principle of conservation of energy. 

Changing v to — ^ and solving equation (3) we get 

ai 



J^sin-lJL + C,. (4) 

If we begin to count time from a position when 3; = 0, we have 



=v- 
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sin— 1 J--= and, consequently, under these conditions C^ = 0, and 



t 

y 



= ^4.sin-X^. ,5) 



—p^t is equal to 0, ^, 2'ir 3^, etc., y = 0, hence 



r' = 



M 



'='V^ <^> 



gives the time of a single vibration, that is the time required for the 
mass to move from the middle of its path to one extremity and 
back again to the middle. The time of a complete period is 



T = 



= 2n^. 



M 



P 

Instead of counting time, from the position where y = 0, we 
might use y = a as the initial position. In equation (4) Ca becomes 

^■^ sin -11. . 

Hence 

. = V4^(si«-l -f _ J|.) = ^I^ eos-l-f . ,8, 

y = a cos-^-^. (9) 



EXAMPLES. 

1. Prove, by differentiation, that equations (8) and (9) are 
solutions of equation (1). 

2. What is the time of a complete vibration of a 2-pound 

mass attached to a helical spring, if a force of 32 poundals stretches 

the spring one foot? 

Ans, 1.57 sec. 



3 D K 



CHAPTER IV 

LINEAR DIFFERENTIAL EQUATIONS OF HIGHER ORDERS 

WITH CONSTANT COEFFICIENTS. 

28. A linear differential equation of the «'* order with constant 
xoefficients is of the form 

-Ai, As, . . . An being constants. 

THE^ OPERATOR D. 

29. In problems of this sort, it is customary to represent the 
t)peration of differentiating and dividing by the increment of the in- 
•dependent variable, by the letter D. That is, when x is the independ- 

«nt variable, — - — is equivalent to D, , „ is equivalent to If etc. 

ax dx^ 

The letter D when used in this way, is called the symbolic operator. 
Using the symbolic notation, equation (1) becomes 

(Dn -i^A^Dn-iJ^A^Dn-^^ . . -\- An)y = /{x). (2) 

An equation such as— r^ + 3-fl + 2y = 4:*, 

dx^ dx 

is represented in this notation by 

{IT + 3D + 2)y = x". 

Expressions written in the symbolic notation may be factored 

like ordinary algebraic quantities. Take first the special case 

(D^ + 3D + 2)y. 

If D were an algebraic quantity, the factors of this would be 

(D + 2)(Z^ + l)y. 

We will show that, when D is the symbolic operator, 

(D + 2) (D + l)y = {D' + 3£) + 2)y. 

■{D-\-\)y is equivalent io-j^-^-y. The expression, (D + o), pre- 

dx 

fixed to any quantity means that the operation -^ is to be performed 

on the quantity and a times the quantity added to the derivative. 
If, in this way, we perform the operations indicated by {D + 2) 

*•" (■£ + y) '* «** 

which proves the proposition for this special case. 

34 



DIFFERENTIAL EQUATIONS. 36 



EXAMPLES. 

Prove by performing the operations 

1. (Z)-2)(D + 2)y=(D^-4)y, 

2. {D —a) {D — h)y = {D^ — {a + h)D -[- ah)y. 

3. (D-a)(D-&)(D-Oy=^-(« + * + 0-S- 

dar ax* 

+ {ah + ac-\-hc)-^ — abcy, 

ax 

Example 3, above, ic a general proof of the proposition for 

expressions of the third order. The terms a, h, and c are the 

roots of the general differential equation of the third order of the 

form of equation (2), regarding D as an algebraic quality, aiid 

omitting fix). 

SOLUTION OF THE DIFFERENTIAL EQUATION WITH CONSTANT 

COEFFICIENTS. 

30. A linear differential equation of the second order with 
constant coefficients may be written 

{D^a){D^h)y = f{x), (1) 

where a and h are the roots of 

Z>^ + ^iD + ^2 = 0, 
regarding D as an algebraic quanttiy. 

In the first member of equation (1), let us represent (Z) — h)y by a 
new variable z. The equation then becomes 

{D^^a)z = f{x) (2) 

This is equivalent to the linear 

_^? a2 = f(x), 

dx '^^' 

the solution of which is 

z = e<^x \e—axf{x)dx + C^e(^x (3) 

Now substitute for z its value, {D — h)y, and get the linear 

From this we get 

y = ^x(e^<^'-b)x( e-<ixf{x)dx^ + Cj^^ + C^e^x ^ (4) 
In a similar way, in an equation of the form of 

everything of the first member to the right of (D — a) is repre- 
sented by a single variable for the first step. Then, for the second 
step, (P — c)y is represented by a second variable. These two 
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operations give a result of the form of equation (4). The final 
solution is 

y = Cy^eax -f C^ebx + C^ecx 4. 

ecxCeib - c)x Ce[a - b)x Ce-ax/(x) dx^. ( 6 ) 

Equations of higher order may be solved in a similar way. 

The last term of equation (5) is called the particular integral 
while those terms which contain arbitrary constants, are together 
called the complementary function. If f(x) is zero, the particular 
integral vanishes, and the solution is the complementary function 
alone. 

THE COMPLEMENTARY FUNCTION. 

31. The method of finding the complementary function is seen 
at once from equation (5) of Art. 30. In an equation of the «'* 
order it consists of n terms of the form 

Ci^ax -j- C^e^^ + Ca^-*" + etc., where a, d, ^, etc., are the roots 
of the algebraic equation (2) of Art. 29. 

EXAMPLES. 

Solve the first two examples by the method of Art. so Tvithout 
using equation (5). 

1. (D' + SD + 2)y = e^, y = C^e-^ -{- C^e-x ^ £1, , 

o 

To verify the solution, perform the operation (i>* + 3Z) + 2) 
on the value of y. 
(rfi + ZD-\-2)y= (4 — 6 + 2)Ci<?-2^ + 

(1 — 3 + 2) C^e-^ + (1 + 3 + 2)^ = e^. 



3. 



4. 



y =. C^^^ + C^e"^ + C^e». 
4^ = a^y, y = C^eax + C^e-ax^ 



CASE WHEN I> + ^iZ2«-l 4. /|jZ?«-2 +....+ ^« = 

HAS BQUAI, ROOTS. 

32. If a pair of equal factors occur we get for the particular 
integral a pair of terms of the form 
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which is equivalent to 

Ci + C2 together form a single arbitrary constant and we have 
then from an equation of the «''» order only {n — 1) constants, 
whereas the number of arbitrary constants should correspond with 
the order of the equation. We will therefore investigate to see 
if equation (5) of Art. 30 applies to this case. 
Take the equation 

(D-a)='(Z)-O3' = 0. . (1) 

For convenience we will arrange the terms so as to bring the equal 
factors (D — o) last. 

(D — 0(^ — a)(^-o)3' = 0. (2) 

If we let (Z? — a)y^=zz, equation (2) becomes 

{P — c){D — a)z^^, (3) 

Since equation (3) does not contain equal roots it may be solved 
by means of (5) of Art. 30 giving 

z^ C^ecx ^ C^eax^ (4) 

ax 
The integrating factor of equation (5) is ^— «■*", hence 

ye-^^ = ^\c^e{c-a)x + c^dx + C^, (6) 

y = C^ecx + \,C^x + C^)eax^ (7) 

We find then that where there are a pair of equal roots the co- 
efficient is of the form 

C\X + C2. 

If we observe equations (5) and (6), we find that multiplying 
by the integrating factor causes the exponential term following C\ 
to vanish, and that this is the reason equation (5) of Art. 30 fails 
for the case of equal roots and for this case only. 

{A similar case is that of the failure of the formula for J x^ dx 
when n = — 1). 

EXAMPLES. 

1. Prove that when three equal factors occur in Z>« -f- A^D^ — '^ 
+ ^jZ?» — 2 -f- . . . -f A^ = the corresponding coefl&cient in 
the complementary function is (Ci;r^ + C^x + Cs). 

Solve and verify results 

2. -g--^+>- = 0, y=(G;r + C.)e*. 

3. (L>' + 2£» + l)y = jr y=(Ci;r + C>)^-* +^ — 2. 
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{Use (4) or (5) Art 30 to And the particular integral. While 
for simplicity we worked out the case of equal roots with f{x) =0, 
the results apply when f(x) does not equal zero. It may he seen 
from Art. 30 that the complementary function does not affect the 
particular integral. 

4. ^--§ffy. + 2y=0, y={C^x+C^)e^ + C^^, 

djT dx 

CASE WHEN I> 4- AiDn — 1 -f- 4^Dn — 2 -^ . ... ^ An = 0, 

HAS A PAIR OF IMAGINARY ROOTS. 

33. Imaginary roots occur in pairs in the form a ± ib. The 

corresponding terms of the complementary function are 

Cie {a^-ib)x _(_ C^e{a-ib)x, (1) 

From analytical trigonometry we have 

gix 4- e—^^ . e^^ — tf— «•«' 
coar = 1 -L- , sin^ir = _ , 

2 2t 

hence 

cosjr + i sinjT = ^«^, cosjr — i %mx = e—*^. 
Cieia + td)x^ Cj^e{<^-ib)x = eax {C^eibx + C^e-ibx) = 

^^ [ ( Q + G ) <^o^x + ( Ci— C, )i sinbx], ( 2 ) 

If we let Ci + C2 = C'i and (d — C2)» = C'2 we get in place of 

the complex expression of equation (1), the simple trigonometric 

form 

eax (C'l cos&.r + C\ sinhx) (3) 

in which C\ and C\ are the arbitrary constants. 

In order that C\ and Cz may be real Ci and Ca must be com- 
plex, so that C\ + Ci is real and d — C2 imaginary. 

As we do not use the original Ci, d, we may drop the primes 

in (3) and write it 

^•^(Ci cosbx + Ci sinbx). 

EXAMPLES. 

Solve and verify the results. 
1. (Z>« — 41) + lS)y = 0, y = ^{Cx cos3;r + C2 sin3^). 



£y_ 
dx" 
Derive the answer of example 2 from that of example 1, Art. 



2. 4^ + &3' = 0, y = Ci cos V <^ ^ +C2 sinV ^ ^. 

djr 



27. 



3. ^— IGy = 0, 3; = Cxe'^ + de-^ + Cs cos2;r + C* sin2jr. 

THE PARTICULAR INTEGRAL. 

34. While the last term of equation (5) Art. 30 shows the form 
of the general expression for the particular integral, the operations 



DIFFERENTIAL EQUATIONS, 39» 

there indicated are frequently difficult to perform. We will there- 
fore give some special methods which greatly reduce the labor. 

THE INVERSE OPERATOR. 

36. If we write Z>« -^ AJ)n-\ •\- A^n--2 ^ AzD^-^ -\- . 
. . . + An as /(D), the linear equation with constant coefficients 

becomes 

f(D)y = X 

where X is any function of x. We will use X in this connection 
instead of f(x) to avoid confusion with f{D), The particular in- 
tegral is then written 

X 

IZdT' 

The expression is called the inverse operator, to distinguish it 

from the direct operator , f(D). The inverse operator indicates 
those integrations and multiplications which are represented by the 

last term of equation (5) Art. 30. For example, -^= J xdx; and 

— is the particular integral of the equation (D — a)y = X. 

U — a 

The particular integral of a linear equation with constant coeffi- 
cients when subjected to the direct operation, f(D)y gives, as we 
know, the function of x of the differential equation. That isi 

represents the particular integral of the equation 
fiD)y = X, and.^^=X. 

OPERATORS COMMUTATIVE. 

36. Direct and inverse operators are commutative. It is imma- 
terial in what order the operations are performed. 

Let (D — a) be a direct operator, and — — =- an inverse oper- 

D — b 

ator. We are to prove that 

From (4) Art. 30, we have 

~—T = ef>^(e-bxXdx. ( 2 ) 

Hence, operating with {D — a), 

{D — a)y^L^ ={b — a)eb^ Ce-^^Xdjc + X. (3) 
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Operating with -=r j^ on (DX — ax), we get 

^n^d^ ^ if *^ Jir-*^Z?^rf^ — aef'^ j e-bxXdx. ( 4 ) 

Integrating by parts, the first term of the second member 

Ce-bxDXdx = e-bxx -f b^e-bxXdx, ( 5 ) 

Substituting in (4) 

^^~f^=^+ {b — a)ebx^e^xXdx, (5) 

which is the same as the second member of (3). 

As any function of D may be made up of factors of the form 
of (D — a), the proof may be extended in this way to functions of 
any degree. 

A special case of the above is that in which a = b. The second 
member (3) reduces to X, thus agreeing with the last equation of 
the preceding article. 

RESOLUTION OF INVERSE OPERATORS INTO PARTIAL FRACTIONS. 

37. Inverse operators may be resolved into partial fractions in 
the same way as algebraic quantities. Consider the expression 

1 

(D — a)(D — b) 
which regarded as an algebraic quantity breaks up into 

a — d\D—a D—bf' 
We are to show that 

1 y ^ 1 / 1 _ 1 Vy 

{D — a){D — d) a — dXD—a D — b) 

operating on (1) with (D — a) we get 

Operating on (2) with (D — b), 

an identity which proves the proposition for the one case. 

EXAMPLES. 

1. Show that-TT^ r-z,X is equivalent to — 



(1) 



(£)_^fl)2 ^ (D — ay D-^a 

2. Show that -JP^ is equal to ^^ , + ^ 



X)=_l '^ ^^ D — 1 D-hl 
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It will help to make our ideas of this subject more concrete, 
if we perform these operations on some simple functions of x, and 
verify the identities. Take example 2 above, and let X = ;r. 

2Dx _ 2 ... 

D^^\ . z)^— r ^ ' 

which is the particular integral of the equation 

(D^-l)y = 2. (2) 

using equation (5), of Art. 80 to find this particular integral, we 

get 

2_ 

In like manner 



_ = 2e^Ce-'^Ce^dx* = — 2, ( 3 ) 



"" +_£--=-^-l + ^-l = ~2, (4) 



D—1 ' £>+! 
whfch is identical with the result of (3). 

EXAMPLES. 

Obtain the particular integral by decomposition of fractions 
and verify the results. 

1. (D^-5D + 6)y = ^, jK=Ci^3^ + C^^+-|: + A. 

THE PARTICULAR INTEGRAL FOR SPECIAL FORMS OF f{x). 

38. It is not always easy to find the particular integral by the 
methods already given. To some forms of X which fortunately 
are those which occur most frequently in physical problems, speciarl 
methods may be applied with advantage. 

CASE WHEN X =^e^^. 

39. Since the nth derivative of e^x^ is a» e^x^ 

f{D)eax =f(a)eax^ (1) 

f{D)eax ^ /{a)eax .ox 

/{£>) fKD) ' 

treating both sides of the equation with the same inverse opera- 
tor. The first member of (2) is equal to ^^, hence the second 
member, 

/('"^^ = eax, 
/{D) 

1 1 ^^ 

and . _ ^^ = 



/{D) f{a) 
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Take the example 



the particular integral is 

Z>« + 4 - 9-h4 l3"' 

and the complete solution is 

> = Ci cos2;r + C2 sin2;r + —=--. 

If it happens that a is a root of f{D) =0, then /(a) =0, and 
= 00, and the method fails. The method of treating a prob- 



f{a) 

lem of this kind is seen from the following example. 

£?f — Q (ly 

(D — 2)(D'— 1) 

Since a = 2 the term {D — 2) renders f(a) =^0. Accordingly we 

will split equation (1) up into two terms 

1 e^ /OV 

{D—2) ' E^— 1 ' ^ ' 

The second term, by the method of this article, gives 

eix __ e*^ 

D^— 1 3~* 

Equation becomes then 

^ 1 ^2;r=£^ (3> 



3 (Z? — 2) 3 

We see then that when (£) — a) is a factor of i(P) it should 
be taken out, the remaining factors treated by the method of this 
article, and the result thus obtained subjected to the operation 

•=r according to the methods of the preceding articles. 

u — a 

EXAMPLES. 

1. (D + 2J9 + 1)3; = ^^; 3;=(C.r + 0^-^ +-^- 

4 

2- 4^- J' = 2. y = Ci^^ + C,^* - 2. 

3. _g.-y=^^, y=(C^J^SLyxJ^C^e-x. 



6 • 12 
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6. (£)'— l)y = cosh2;r. 



y = C^e^ + C^e-x 4. 



1 ^ J^ e-'ix 



3 



ory = Cj^^ + C,^-* + -i- cosh 2Jir. 

7. (D'- l)y = ^^, JK = C^e' + C,e-x - €1 . 

8. (D — D'y = cosx, y = iCxX+C)e* - ^. . 

l^.Make use of the relation cos;r= -£ — it Y 

9. -^ + 4-^ + 3>^ = cos 2x, 

y = C^e-x + Q^S;^ + 8 sin 2;r -. cos 2^ ^ 

65 

CASE WHERE X = slno^T or cosa^r. 

40. D sxnax = a cosax 

£)* sino.r^: — a* sinajr 
Hence f(D') sinax = f( — o*) sinojr. (1) 

Operating on equation (1) with we get 



f(D') sinax _ f( — a^) sinax 



Since the first member of (2) is equal to sina;r we have 



(2) 



f(^a^) sino^^^j^^^^ (3^ 



Dividing (3) by the algebraic quantity f( — a*) we get 

singjr __ singer 

Similarly 

coso;r _ cosajt 

Take as an example 

(D* — 4)3; = sin;r 
the particular integral of which is 

sin;r __ sinjr sinjg 

/>* — 4 —1 — 4 5 ' 

The complete integral is 

y = Ci^^ + Cie-^^ — 



sin^ 



6 

When, as is generally the case, contains factors of the first 

1 V**"^ y 
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degree of the form -j^ these may be brought to the second 

degree by operating on numerator and denominator with (Z) + 4r). 
For example 

sin2.r __ {D + l) sin2;r _. 

{D" + 3) {D —I) (D= + 3) (D^ — 1) 

(D + 1) sin2;r _ 2 cos2;r + sin24r 
5 6 

If quadratic factors occur in containing both D' and D. 

D* may be replaced by — a^ and combined with the constant term 
and the resulting linear factor treated as in the example above. 
For example 

cos2;r __ cos2.r c os2;ir _ 

D^ + 2D + S — 4 + 2D + 3 2D — 1 

(2D-\-l) cos2.r _ 4 sin2jr — cos2;r 
W — 1 17 * 

That this is the correct solution is seen by performing the 
operation (Z)^ + 2D + 3) upon the final result. 

When ^^ ^^ contains the factor -— ;— — ^ the method fails for 
f{D) D' + a' 

= 00. To evaluate — — j — ^ we will give to a in sinajir 



-^a' + a D- + a 

2l small increment h. The expression becomes then 

sin (ax -j- hx) _ sin(ajr + h,r] _ 

D' + a' — (a' + 2a}i + h') + a' 

sin ax cos hx + cos ax sin hx / •• % 

2a/! + h' * ' 

As h approaches zero as a limit, cos hx becomes practically unity, 
sinhx becomes equal to hx, and h^ is negligible. 

We get then from ( 1 ) - ^^ ±cosax__ ^2) 

2a li 2a 

The complementary function of 

(D^ -i-a^)y = sina;r 
is Ci cosojr + Ci sinax. The second term of this complementary 
function has the same variable factor as the first term of the 
particular integral. Since the constant Cs may have any value from 
+ 00 to — 00 they may be combined into a single term. The com- 
plete integral of 

(D' -{- a^)y =^ sinax is then 



Ci cosa.*" + Ci sinax — 



X cosax 
2a 



as may easily be proven by performing the operation (D^ + a*). 



DIFFERENTIAL EQUATIONS, 4^ 

In a similar way it may be shown that 

cosfl.r _ X sinojr 

EXAMPLES. 

Solve and verify results. 

1. -^-9y = cos^, ^ = Ci^3a^ + C,^3a^ - -£?H. . 

ax* 10 

dx^ ax 

y = C^e'^x + c^e^^ + cos^ + sin^ ^ 

3. 4^ + 4y = sin2.r, y=^C^ cos24r + G sin2.r ^JL^^^L . 
ax^ 4 

3^ = Ci cos2.r + C, sin2;tr + C^e-x + 3cos3x - sin3^ 

5. (Z>2 __ 7£) -|_ 12)3, = e'^ + sin3.r , 

1/ = C^e^x + (r,^4a? + £?? 4- sinS^r + 7 cos3;r. 
* ^ * ■ 2 ^ 150 

6. *1 + «y = cos^^, y = a-aa, + ^ ^^"^^ + ^ ^Q^^-^ 
^;r &^ + a* 

7. Integrate ^«^ zo^hxdx by means of inverse operators. 
The particular integral of the linear of the first order 

-Z + ay — cos&jr (1) 

dx 

is expressed by e—iif> j eo^cosbxdx. 

If we equate this to the expression obtained with the inverse 
operator we get 

(T-M C eaof cosbxdx = -|os&£^ ^2) 

^ D -{- a 

Multiplying (2) by e<uf 

Ceaxcosbxdx = ^aa,.£0s6^ ^ ^ ^^ b smbx + acosbx^ 
•^ D-\-a D +a' 

omitting the constant. 

8. Solve example 9 Art. 39. 

n ci,« , *u»4. coshajr! _ cosha^r _„ j sinha^r __ sinhojr . 

9. Show that -^^^^■-^^^. and -^^^^--^^^ 

cosh2^ 



10. (D" — Dy = cosh2;ir, y = C^e^ + C^e-^ + 



3 
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CASE WHEN X = X^y WHERE m IS A POSITIVE INTERGER. 

41. To find the value of , ^ , expand — — — into a series in 

ascending powers of D^ and operate with this series on jc^. Since 
w is a positive interger, D^x^ is a constant, and all terms beyond 
this are zero, so it is only necessary to carry the series to /?*». For 

x^ 
example evaluate 



U' — l 
1 _ 1 



= — 1 — D» — £)*, etc. 



(1 + D' + D')x' = — (x^ + 6 x+0) = — jr» — 6;r. 



EXAMPLES. 



1. ^2 + ^->' = ^' 



y = Ci cos 2.r + C2 sin 2;r + -L^ ^r* — 5.r»^-i^^ . 

2. (D' + SD + 2)y = 2 + x\ 

3. ^-\.ay = x^, 
ax 

y ^ ce^.^ x^_4x^,12x^_24^,24^^ 

4. Integrate e^^x^dx by means of inverse operators. 
As in ex. 7, Art. 40, we find that the integral is 

D ^a ~ \a a« "^ a* / * 

5. (D» + D)y = x% 

y = Ci^-^+G+-^(:r* — 4;r» + 12Ar'^ — 24.r + 24^ . 

D"' + D D(l + D) Z>V / 

We may complete the solution in several difierent ways. We might 

perform the operation — -. first and then the operation 

(1 — D + D^— . . ). 

J_;r3 = ^, and 
V 4 ' 

(^l — D + D' — D' + D* . . )— = 

_L^;i'4 _ 4^3 ^ ^2x' — 24.r + 24^ . 
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Again we might combine the factors into (—- — 1+^ — ^D*+^ • • ) 

and operate with this on ;r*. Finally we might perform the opera- 
tion (1 — Z> + £>* — £>*+ . . ) first and then integrate the re- 
sult. 

CASE WHERE X z=. e^^ V , V BEING ANY FUNCTION OF X. 

42. Take first the case where f{D) is of the, first degree. 
^ is the particular ^'ntegral of the equation 



{D--a)y = ef^^y. (1) 

Assume that y=^e^^s, and equation (1) then becomes 

(D — a)c^^z = e ^^ V, (2) 

Performing the operation (Z> — a) on e^^z in the first number of (2), 

enix(p — a + m)z = ^ w^ F. (3) 

Dividing by e^^ ^ and performing the operation -=r ; 

D — a-|- w 

(4) 



D — a-\-m 

We may now multiply by the factor e^^ and substitute for e^^z 
its value y. This gives 

y = ^^-n — ^ — (5) 

D — a-\- m 
When f(D) is of the second degree, we have the form 

(D — a) (D — b)y = e^^V. (6) 

We may regard (D — b)y as a single variable u, and write 

(D — a)u = e^^V. (7) 

Applying the results of equation (5) to this 

u = e^^— — ^1— = (D — &)y . (8) 

D — a-f- m 

Regarding •=; j — as a single function of x and applying equa- 

V — a -\- tn 

tion (5) again we get 

^ iD — b + m)(D — a + m)' ^^ 

We proceed in a similar way when f(D) is of higher degree. In 
any case we find 

f!!^ = e^-' ^ . (10) 

EXAMPLES. 

1. {£^ + 2D + l)y = e'i^cos X. 

y = ( C^x + Cg ) e-^ + ^{ 4 cos ^ + 3 sin ^ V 
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2. ( Z?2 — 4Z? + S)^ = ^2;rcos X, 



y =: de^^ + de^ — 



3. £y--y^e-'^, , 

y = C.e^+C.e-^~£:^(^x' + 2x' + ^x' + ^x^-^^ . 

4. Integrate e^^ cos Sxdj:: by means of inverse operators. 

5^ 5 cos 3.r + 3 sin Sx 
34 

5. J^2;r;i:«^;,;, ^(^Ax^ — 6x^ + 6x — S\ . 

6. -:^ — y = X cosh ;r, 



JV 



= Ce^ -\- —( x^ex — 1xe-x — <?-:r \ 



7. -7^ = .r^ sin ;r, y^^C — :^ cos ;jr + 2^ sin ^ + 2 cos x, 
ax 

8. (Z>2_4)3, = ^ cos ^, 

y = Cie2^ + ^8^-2^—^ COS 4: + ^ ;r sin ;ir ■— j^cos x, . 



PROBLEMS. 

1. Solve the problem, Art. 27, by the symbolic method. 
From the equation 

M^ =^Pyvfc get 
or 

y = Ci cos ^— / + Ca sin-y/— A (1) 

When ^ =0, y = 0, hence Ci = 0, and 

y = C2 sin ^-^^» (2) 

When sin \-n^^ is unity the value of y, which is then the max- 
imum, is equal to Cz, Writing a for this maximum y we have 

y=asin^-^, (3> 

2. Euler's formula for the strength of a column. 

Fig. 5 represents a column or compression member (drawn hor- 
izontal for cenvenience) subjected to a compression at the ends. 
The moment at any section distant x from the left end is — Py, 
where P is the pressure at each end, and y is the displacement of 
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Figure 5. 
the neutral axis of the section. When y is positive -^ decreases 

as X increases, and — ^ is negative ; and when the deflection is down- 

dx 

ward making y negative, -j^ increases as x increases, and — ia 

positive. Hence in either case 

EI^=-Py. (1) 

dx' 

y = Ci cos -v/-£7^+ ^« sin aZ-^v*^* (2) 

Putting in the limiting conditions wc fmd that y = ^ sin aZ-etv*" » 
where A is the maximum deflection, and th2it^-=r-l = rnr, n being 

an integer. From this we get Euler's formula, P = — j^ — 

When the column bends in a single curve, as in the figure, n is unity. 
Column with eccentric load. 
When as in Fig. 6, the load is not applied at the center of the 



Figure 6. 

cnJ-. of a column but at a distance k from the center, the differential 
equation is 

E/^.= -P(k+y}, (1) 

IP IP 

y = Ci cos^-^jix: -\-C2 sin ^-^x — k., (3) 

From the condition that 3^ = when x = we fir. d Ci = k. Hence 

y = k cosyj-j^x + CiSm^—x — k. (4) 

4 D E 
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From symmetry it is evident that the maximum deflection is 

sl^ the middle ; hence -7^ = when x = — -.. 

dx 2 

J)ifferentiating (4) and substituting we get 

Substituting in (4), combining, and transforming, 

y = *[sec^^4cos^5(^-4-)-l]. (6) 

The maximum deflection at the middle where x = —— is 

2 



A = *[sec^'£4-l]. 



^See Merriman, Mechanics of Materials, p. 132.) 
Beam -and one end fixed and with a force applied in any direc- 
tion to the other end. 

Fig. 7 shows a beam with the right end fixed horizontally, and 



Figure 7. 

"with a force P, making an angle a to the left of the vertical, ap- 
plied at the left end. 

We take the origin at O, the position of the end of the deflected 
beam. Resolving the force P into its horizontal and vertical com- 
ponents and multiplying these by y and x respectively, we get the 
.<liflFerential equation 

EI^ = Py sin a — Px cos a (1) 

dx^ 

[^^_^], = _^^.cosa (2) 

\P sina \P sina 

y =zC^e + C^e + ^ cot a, ( 3 ) 

"Using the conditions that j' = when ;r = 0, and —^ = 0, when 

ax 

jr = /, we get 
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Ca — — w — 



cota 



V 



^ ^^"^ cosh ^^iJ2^/ ' *** 



EI \ EI 

Substituting in (3) 

•„!, /P sina 

y — zo\.o.m x^ ^ ^' 



VP smo. ^^„i, IP sina. 



{See Church, Mechanics of E^ngineering, p. 357). 

If the direction of the force P is to the right of the vertical 
making an angle i3, the moment of the horizontal component is re- 
versed and the differential equation becomes 

Putting in the limiting conditions we get 

sin '^^ 

^^ -X 1. (3) 



y = cot/3 I 



W^^cos^f^/ 



In these problems we have the right end of the beam horizontaL 
The formula will apply to a beam in any position if the X axis is 
taken parallel to the fixed end, and the angles a or P be the inclin- 
ation of the force to the corresponding Y axis. 

ELECTRICAL PROBLEMS. 

Circuits having resistance, inductance and capacity. 
45. The E. M. F. equation is 

^i + L§ + ^ = e. (1) 

the positive direction of the current being towards the condenser. 
Substituting and multiplying, 

^^3 + ^C^ + , = Ce. (2> 

Whatever the value of e. the complementary function is the same^ 

The roots of LCD' + RCD + 1 = being ~ ^^ -^ ^'^ ~ ^^^, 
the complementary function is 

— RC+ y^2C2— ALC —RC—VR'^C^— 4LC _, 



C^ l"-^ + C,e ^^ , (3» 

when V IVa — ALC is real. 
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When y R^C^ — 4LC is imaginary, the complementary function 
is 

e-2L[C,^^^ ^^ 1 + C^cos^^ ^^ tj (4) 

When R^O = iLC, we have the case of equal roots and the com- 
plementary function becomes 

{Ci -f- Cj/)*? 2L. (5) 

The particular integral in any case is 

Ce 



LCD' + RCD + 1 



(6) 



Discharge of a Condenser. 

In this case e = and the complementary function is the com- 
plete integral. We have then only to find the constants for each 
of the cases above. 

Case of real, unequal roots: /?*C*>4LC 

_ /gC - |/y?2C2 - jLC _ ^C -f VR2Ci — ALC 

^ ^LC ^ 2LC 

We may leave the equation in this form, or we may transform it 
before determining the constants. In the latter case 

, = .-^(c',sinhil!^^C, + c'.coshJ^^i?£^)(2) 

where C\ = ^'7^' and C'2 = ^^ + ^' • 

From the conditions that when * = 0, q=Q, and « = we get 

C,^Q and C,= ^ ^,^_^^^ . 

q = Qe-^( ^^ sinh V^'C'-4LC ^ 

^ ^ \iJR^C — iLC 'il.C 



Equation (3) may be reduced in the same way as a similar 
trigonometric function to 



ftanh-tV^!C^=«£.),4)- 



To get the current we differentiate one of the expressions for 
q, preferably (3). The result reduces to 
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t = — — === e 2L Sinn ^ /, ( 6 ) 

Ca^^ 0/ imaginary roots : R^O <iiLC, 

_J?i/^ . V4LC — R^'C . , /^ V4LC — R^C . \ ,a^ 

q - e 2Z^CiSin ^ZC "^ Qcos ^^^^^^^^ / ^, (6) 

Show that C2 = Qi* 



^\LC — R'O ' 
. ( 



?=-???^=^i^^ 2^ sml 577= ^ 



V4LC-je'C' V 2y:c 



7) 






2^ .-^siuJ^HgP!^/, ..8) 



y?c / 



V4LC — /e^C 2Z.C 

Since (7) and (8) are sine functions i and ^ may have, 
positive and negative values. The discharge is said to be oscilla- 
tory. In the preceding case, since the hyperbolic sine is always 
positive, the discharge is seen to be non-oscillatory. 
Case of equal roots : R^O = 4LC. 
The integral is 

q=z(Ci + C2t)e 2^. 
Show that Ci = Q, - 

jRt 

IT /. , Ri 



,=Qe -(!+-§-> 



< = — tf 



Rt 

2L QR^t 

4Z,« 



and since Q = CE, and -4— = 1, 

4Z* 



Rt 



L 

Charge of a condenser with a constant E. M. F. 
The charge of a condenser with a constant voltage E differs 
from the discharge in only two respects: at the initial time 
g = 0, and there exists the particular integral 

CE Q _^ 

LCD^ + KCD + 1 LCD' + RCD + l ^' 
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Show that for the case of imaginary roots 

^ ^\ V4LC— ie*c? \ 2Z<r 



RC J ' 



Derive the expression for q for each of the other two cases. 
Derive the expression for i for each case. In what respect does 
the current of a charging condenser differ from that of a discharging 
condenser? 

Circuits in which e = £ sintot. 
The particular integral is 

CE sinttf/ _ CE sinto/ 

LCD" + RCD + 1 RCD — (LCo^' — 1 ) 

_ CE(RCD + (LC<o' — 1 ) ) sin«/ 



R^CD' — (LCo^^ — 1) ^ 

_ Rcos(ot-{- FLw — —-I sinw^ 

_ jD^ L CtaJ 



(4) 



A short time after the circuit is closed the terms of the comple- 
mentary functions vanish and (5) gives the value of q. The cur- 
rent after a short time is expressed by 

,• _ ——===== sin I a,*-tan -il"^ - 74-11' <»> 



\ 
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